Tail states in clean superconductors with magnetic impurities 
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We analyse the behavior of the density of states in a singlet s-wave superconductor with weak 
magnetic impurities in the clean limit. By using the method of optimal fluctuation and treating 
the order parameter self-consistently we show that the density of states is finite everywhere in the 
superconducting gap, and that it varies as In iV(.E) oc — \E — Ao|' 7-< ^' 4 near the mean field gap edge 
Ao in a d-dimensional superconductor. In contrast to most studied cases the optimal fluctuation is 
strongly anisotropic. 
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One of the most intriguing problems in modern con- 
densed matter physics is the combined effect of correla- 
tions and disorder on the ground state and the electronic 
properties of solids. Interplay of superconductivity and 
impurity scattering is an example of a situation where 
such an effect is relevant experimentally, and where it has 
been studied for over fortyyears, beginning with the sem- 
inal papers by Anderson |lj and Abrikosov and Gor'kov 
(AG) |2(. Despite such a long history, however, the spec- 
tral properties even in the simplest case of a singlet s- 
wave superconductor continue to attract attention. 

Since nonmagnetic impurities do not break the time- 
reversal symmetry, the "hard" gap, Ao, in the single 
particle excitation spectrum is unaffected by weak po- 
tential scattering [lj. At the same time scattering by 
magnetic impurities destroys the phase coherence in the 
superconducting state, and, consequently, leads to the 
suppression of the gap, and of the superconducting or- 
der parameter. AG analysed magnetic scattering using 
self-consistent Born approximation, and concluded that a 
hard gap in the energy spectrum survives up to a critical 
(average) concentration of weak magnetic impurities; it 
is followed, upon increasing impurity concentration, by a 
narrow region of gapless superconductivity, and then by 
destruction of the superconducting condensate. 

After some early work 0] , resurgence of interest in this 
problem started with the paper of Balatsky and Trugman 
[4|, who argued that rare regions, where impurity concen- 
tration is sufficient to locally destroy superconductivity, 
yield a finite density of states (DOS) at the Fermi level. 
Their work was followed by other analyses Here we 
elucidate the nature of the subgap states, and provide a 
quantitative analysis of the energy profile of the DOS in 
s-wave superconductors with weak magnetic impurities. 

In the AG theory the effect of disorder is controlled by 
a dimensionless parameter, At s , where t s is the scatter- 
ing time due to magnetic impurities, and A is the value of 
the superconducting order parameter. In particular, the 
single particle spectral gap is A = A(l — (At s )~ 2 / 3 ) 3 / 2 , 
indicating the onset of the gapless superconductivity at 
At s = 1. In the regime studied here, Ar s ^> 1, the AG 



theory predicts a gapped quasiparticle spectrum. 

These results are obtained by carrying out a standard 
impurity averaging procedure. It is clear, however, that, 
among all the realizations of the impurity distribution, 
there exist regions where the resulting potential generates 
localized quasiparticle states at an arbitrary energy be- 
low the gap edge. Such localized states were extensively 
studied in doped semiconductors 0, 0] • For a particular 
energy, E, the most probable (albeit still rare) configu- 
ration of impurities that creates a state at E, and there- 
fore contributes the most to the DOS, N(E), is called the 
optimal fluctuation (OF) 0- Such rare regions provide 
nonperturbative corrections to the DOS determined in 
the framework of self-consistent Born approximation. 

We employ the OF method in a singlet s-wave super- 
conductor with magnetic impurities. We also consider 
the self-consistent suppression of the superconducting or- 
der parameter; it is known that in d-wave superconduc- 
tors it significantly affects the low-energy DOS ||. We 
find that the density of states is finite everywhere in the 
superconducting gap. Just below the AG gap edge, the 
DOS is N{E)/N„ oc exp (-\E - Aof-^^), where N 
is the normal state DOS, and d is the number of spatial 
dimensions. In contrast to other known cases, the OF 
is anisotropic, with its transverse size much smaller than 
the longitudinal extent. 

We consider a singlet s-wave superconductor. In the 4- 

space of the wave functions (ip*(r), ^*(r), 4>]{r), , 
the mean field hamiltonian is 



H = £r 3 + A(r) Tl a 2 + U 



(1) 



Here £ = — V 2 /(2m) — /x is the kinetic energy of a quasi- 
particle with respect to the Fermi level, /i, and Tj and <7j 
are the Pauli matrices in the particlc-holc and the spin 
space respectively, so that Tj<7j is a 4x4 direct product. 

Potential due to impurities, U, includes both potential 
and spin- flip scattering processes. When the potential 
scattering is dominant (motion of quasiparticles in the 
OF is diffusive) properties of the low-energy states were 
explored in Ref. ||. However, there exists experimental 
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evidence that in some situations the magnetic scattering 
is dominant: upon increasing the impurity concentration 
the increase in residual resistivity ratio correlates with 
the suppression of the superconducting transition tem- 
perature 0- Guided by this insight, we consider only 
the magnetic scattering, and expect our results to remain 
valid for as long as it is stronger than, or is of the order of, 
potential scattering. Hence we write U = U(r) • s, where 
s is the electron spin operator, U(r) = J2i JSiS(r—ri), J 
is the exchange constant, and Si is the localized impurity 
spin at a site i. 

The main physical difference between our analysis and 
that of Ref. |5| lies in this choice of the scattering po- 
tential. In the regime studied here the mean free path 
significantly exceeds the coherence length, and hence the 
motion of the quasiparticles in the optimal fluctuation is 
ballistic, leading to a substantially different physical pic- 
ture of the tail states and the optimal fluctuation, and to 
a different energy dependence of the density of states. 

The states with energy E < A exist in rare regions 
where the amplitude of the impurity potential differs sig- 
nificantly from its typical value. Therefore in determin- 
ing N(E) it is sufficient to consider only such configura- 
tions of the impurity potential, for which E is the lowest 
quantum mechanical energy level; fluctuations where E 
is the energy of a higher bound state are exponentially 
less probable 0. Also, in essentially all the energy range 
below the gap the size of the optimal fluctuation is sig- 
nificantly greater than the distance between impurities, 
so that the exact impurity potential can be replaced by a 
smooth function, averaged over regions containing many 
impurities, but smaller than the characteristic size of the 
wave function in the optimal fluctuation [2j. 

Hence we consider an uncorrelated potential with a 
gaussian probability density 

P[U] cx exp (--jjs J ef*rU 2 (r)) , (2) 

where U 2 — n imp J 2 S(S + l)/3 is related to the scatter- 
ing time via t^ 1 — 2-kNqUq, and rii mp is the average 
impurity concentration. We ignore interactions between 
the magnetic impurities: it was shown in Ref. |lDl | that 
the RKKY interaction and glassy behavior of impurity 
spins modify the AG results very weakly. We also do not 
include quantum dynamics of the impurity spins, and 
therefore cannot account for the Kondo effect. This is 
justified either when the Kondo temperature of individ- 
ual impurity sites is much smaller than the superconduct- 
ing transition temperature, Tr- -C T c (and depletion of 
states at the Fermi level prevents screening of the local 
moment), or in the opposite limit, Tk S> T c , when the 
moments are quenched already in the normal state 

The density of localized states in the fluctuation region 
of the spectrum is then 

N(E)= [v\JP[U]6(E-£[U]), (3) 



where £ [U] is the lowest energy eigenstate in the realiza- 
tion U of the impurity potential. For rare configurations 
the integral is evaluated using saddle point approxima- 
tion to give In [N(E)/N ] ss —S[U op t], where the optimal 
fluctuation is obtained by minimizing the functional 

S[V] = 2fP I + X ( £ W - (4) 

with respect to both U and the Lagrange multiplier A. 
The difficulty in minimizing the action is in the nonlinear 
nature of the equations: optimal potential U depends on 
the wave function of the particle in this potential. 

The method of optimal fluctuation allows for a simple 
physical analysis. Consider first a semiconductor. In a 
potential well of depth U (all energies are measured from 
the band edge) and size L the energy of the localized 
state is of the order of U + l/(mL 2 ) = E (h — 1). In the 
optimal fluctuation E ~ U ~ L~ 2 , so that the action for 
such fluctuation is S[U] « L d U 2 /U$, or ln[N(E)/N ] ~ 
-\E\ 2 - d / 2 /U$. This is exactly the result obtained by 
Lifshits, and is confirmed by the solution of the nonlinear 
equations for minimization of action in Eq.Q) fj, 0. 

The difference between a potential well in a doped 
semiconductor and in a superconductor is twofold. First, 
because of particle-hole mixing the hamiltonian Eq. Q is 
a matrix in particle-hole and spin space. Second, we are 
concerned with quasiparticles close to the Fermi energy. 

We assume that ferromagnetic fluctuation maximizes 
the effect of the impurity potential ^3] • Consequently, we 
consider such a fluctuation, and choose the direction of 
the impurity spins along the y-axis, so that U = U(r)<j2- 
Performing a spin rotation, 02 — > 03 in Eq. , we obtain 
a hamiltonian diagonal in the spin space, 

H ± = fr 3 ± A n ± f/(r). (5) 

It is therefore sufficient to consider only one spin orienta- 
tion. Let us again consider the problem qualitatively, and 
concentrate first on the one dimensional case ignoring the 
suppression of the order parameter. We linearize the ki- 
netic energy near the Fermi surface, £ ~ivp(d/dx), 
so that typical kinetic energy in an OF of size L is 
£ ~ vp/L. Then the energy of a quasiparticle in the 
optimal fluctuation (measured from the Fermi energy) is 
E ~ U + x/Aq + Vp/L 2 . For the energies close to the 
superconducting gap, (A — E)/A <C 1, the OF is large 
(L > £ = wf/A ) and shallow (|f7|/A < 1), so that 
E — Ao ph U + Vp/(AoL 2 ). Introducing the dimensionless 
energy e = E/ Ao, we obtain, in analogy with the argu- 
ments above, |Z7|/Ao — £,q/L 2 ~ 1 — e. Notice that the 
size of the fluctuation is indeed L ~ Co/Vl — e ^ Co- As 
a result, we find S[U] « LU 2 /U§ = Agf (l - e) 3/ 7^ 2 - 
From the definition of Uq it follows that 

- In « S[U opt ] ~ (A r s )(l - ef 2 . (6) 
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The energy dependence in Eq. © is identical to the result 
of Lifshits in d — 1 . This follows from the expansion in 
£ <C A : even though £ oc the expansion is in £ 2 . 

We now verify these estimates by a considering the 
energy of a particle in hamiltonian Eq.JSJ) for spin "up" 



£+[U] = (H+) = (<F|£r 3 + A ri + U\9), 



(7) 



where ^Sf is the normalized wave function of the particle. 
Minimization of Eq. (@J with respect to U gives 



U(x) = -\U*(V\ 



~W 



l*>, 



(8) 



while minimization with respect to A dictates -ff+l'F) = 
E\fy). We first ignore the self-consistent suppression of 
the gap, which means U(x) = — XUfi (^* (x)^ (x)) , where 
(\f r *\I f ) denotes the scalar product in particle-hole space. 
Then the Schrodinger equation takes the form 



d 

' iVF dx T3 + A ° Tl ~ AC/ oO**) 



* = EW. (9) 




FIG. 1: Bottom panel: comparison of the optimal poten- 
tial for A = Ao from Eq. dlUt with that obtained from the 
self-consistent numerical solution, for different values of the 
bound state energy E = eAo. Top panel: the corresponding 
self-consistently determined profile of the order parameter. 
Maximum in A(x) at x = corresponds to small gap sup- 
pression at an extremum of the potential, see Eg. (1121 , 



This equation is solved by introducing the bilinear forms 
(x)Ti^ (x) , and yields the optimal fluctuation 



U(x) 
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2A e + cosh^yT^ - ? /£„) ' 
which corresponds to the value of the action 



S[U] = 8tt(A t s ) yl — e 2 — earccose 



(10) 



(11) 



We immediately notice that for e : 



1 the length scale of 
, its depth is U ~ 
(87r/3)(A r s )(l - 



the optimal fluctuation is £o / yl — C- 
Ao(l — e 2 ), and the action S[U op t] — 
e 2 ) 3 / 2 , in complete agreement with our estimates above. 

We now show that the self-consistent suppression of 
the order parameter does not appreciably change our re- 
sult. Self-consistency is achieved by including the varia- 
tion of the gap into the variational derivative 8H + /SU. 
We notice that, at T = 0, uniform U does not suppress 
superconductivity. Consequently, A depends on the gra- 
dient of the potential. For 1 — e -c 1, the potential varies 
smoothly, so that dU/dx <§C Ao/£o can be accounted for 
perturbatively. The leading local correction to the gap is 



ijf. 

6A 2 



dU 
dx 



(12) 



and the correction to the action for e ~ I is 5S ~ 
(1287t/315)(A t s )(1 - e 2 ) 7 / 2 . At lower energies the self- 
consistent A(x) has to be computed numerically. The 
results are presented in Figs. I and 2. ft is clear that the 
suppression of the gap, even for the states with E <C Ao 
is incomplete; since \U\ < 2Aq, the gap remains at a 



significant fraction of Ao throughout the OF. This jus- 
tifies the expansion in the bare energy, £, in our quali- 
tative analysis. Consequently, the optimal action com- 
puted self-consistently differs at most by 10% from that 
computed assuming a uniform gap, see Fig. 2. 

Having demonstrated that the qualitative considera- 
tion are in excellent agreement with the full solution of 
the problem in d = 1, we discuss the multidimensional 
case. In a doped semiconductor the OF in any d is spher- 
ically symmetric 0, ■ This results from the balance be- 
tween lowering the particle energy in a large and deep 
fluctuation, and the probability cost of such an OF. 

Since electrons in a superconductor move with the 
Fermi velocity, the wave function of the subgap state is 
concentrated along the quasiclassical trajectory, which is 
a chord in a potential of any shape. Consequently, there 
is little energy cost in reducing the size of the OF in the 
"transverse" direction, while the smaller volume makes 
such fluctuations more probable. As a result, the op- 
timal fluctuation is anisotropic, and strongly elongated 
in one direction. Choosing this direction as the x-axis, 
we can write the wave function of the subgap state as 
*(x,y) = exp(ikpx)$(x, y), where y denotes the trans- 
verse d—1 coordinates, and $ is a slowly varying function. 
Therefore the kinetic energy of the quasiparticle is 



Akpx 



ox 




The transverse size of the fluctuation can therefore be 
reduced until the second term becomes comparable to 
the first, i.e. L y ~ (XpLx) 1 / 2 , where \p ~ hp 1 is the 
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FIG. 2: Optimal action, normalized by So = 87r(Aor 3 ). Lo- 
cal approximation: A(x) = Ao + SA(x), with SA(x) from 
Eg.ltHl. 

Fermi wavelength. Consequently, we find |J7|/Ao ~ 1 — e 
and L x ~ £o/\/l — £ j and 

5[C/ opt ] ~ A^" 1 ^ * (Aor s ) (^J (1 - e)T , 

(14) 

where £?p is the Fermi energy. Ea. (|14|) is the main result 
presented here. The action for the anisotropic fluctuation 
is smaller than that for an isotropic OF, by a factor of 
(^ jF /Ao) (d - 1)/2 (l-e)-( d - 1 )/ 4 , so that the corresponding- 
DOS is exponentially higher. 

There are two limitations on the validity of the results 
obtained here. First, since the optimal fluctuation is a 
result of a saddle point approximation for the functional 
integral, Eq.|j2J|, it is only valid when S[U op t] ^> 1, or 

l- e »(A r s )^ (Jij 7 " . (15) 

For d — 1 this condition becomes 1 — e ^> (Aqt s )~ 2 / 3 , 
while for d = 3 it does not depend on the gap amplitude, 
l-6> (kpl)^ 1 , where I = vpr s is the mean free path. 
In d = 1 the region of validity is extended by almost an 
order of magnitude in comparison to this estimate as the 
action has a large numerical factor ss 24. 

Second, when the characteristic size of the OF L > I, 
our assumption about the ballistic motion in the fluctu- 
ation is invalid, and a crossover to the diffusive regime 
studied in Ref. occurs for 1 — e < (Aqt s )~ 2 in any di- 
mension. When d = 1 for (Aqt s ) > 5 the saddle point ap- 
proximation becomes invalid before the diffusive regime 
is reached. For d > 1 and typical A /Ep — 10~ 3 , the OF 
method works up to the crossover. Taking (Aqt s ) ~ 10, 
we find that our results hold to within 1% of Ao, while 
the expansion in £ is quantitatively valid for e > 0.9, and 



qualitatively for e > 0.75, providing a significant window 
of applicability for our DOS. 

Experimental verification of our results, and the un- 
derlying physical picture of scattering on randomly dis- 
tributed impurities, requires averaging the tunneling con- 
ductance over regions containing many impurities, and is 
best done at energies just below A. We suggest averaging 
the tunneling spectra, obtained from Scanning Tunneling 
Spectroscopy, over several distinct areas of the sample, 
each containing a large number of impurities. 

To summarize, we have analysed the density of subgap 
states in an s-wave superconductor with weak magnetic 
impurities using the method of the optimal fluctuation. 
We concentrated on the clean limit, ! > ^o, when the mo- 
tion of particles in the optimal potential is ballistic. We 
find that the optimal fluctuation in this case is strongly 
anisotropic, and that the density of states varies as a 
stretched exponential below the gap edge, with a power 
that depends on the dimension. 

This research was supported by the DOE under con- 
tract W-7405-ENG-36, by the NSF Grant PHY-94- 
07194, and by Pappalardo Fellowship. We are grateful 
to the ITP Santa Barbara for hospitality and support. 



[1] P. W. Anderson, J. Phys. Chem. Solids 11, 26 (1959) 
[2] A. A. Abrikosov and L. P. Gor'kov, Sov. Phys. JETP 12, 
1243 (1961). 

[3] I. O. Kulik and O. Yu. Itskovich, Sov. Phys. JETP 28, 
102 (1969); A. I. Larkin and Yu. N. Ovchinnikov, ibid. 
34, 1144 (1971). 

[4] A. V. Balatsky and S. A. Trugman, Phys. Rev. Lett. 79, 
3767 (1997) 

[5] A. Lamacraft and B. D. Simons, Phys. Rev. Lett. 85, 
4783 (2000); Phys. Rev. B 64, 014514 (2001); I. S. Be- 
loborodov, B. N. Narozhny, and I. L. Aleiner, Phys. Rev. 
Lett. 85, 816 (2000). 

[6] I. M. Lifshitz, Sov. Phys. Usp. 7, 549 (1965); Sov. Phys. 
JETP 26, 462 (1968); J. Zittarz and J. S. Langer, Phys. 
Rev. 148, 741 (1966); B. I. Halperin and M. Lax, ibid. 
722. 

[7] I. M. Lifshits, S. A. Gredeskul, and L. A. Pastur, Intro- 
duction to the theory of disordered systems, (John Wiley 
& Sons, New York, 1988). 

[8] W. A. Atkinson, P. J. Hirschfeld, and A. H. MacDonald, 
Phys. Rev. Lett. 85, 3922 (2000). 

[9] M. A. Woolf and F. Reif, Phys. Rev. 137, A557 (1965); 
A. S. Edelstein, Phys. Rev. Lett. 19, 1184 (1967). 
[10] A. I. Larkin, V. I. Mel'nikov, and D. E. Khmelnit- 
skii, Sov. Phys. JETP 33, 458 (1971); V. Galitski and 
A. I. Larkin, cond-mat/0204189 
[11] E. Miiller-Hartmann and J. Zittartz, Phys. Rev. Lett. 26, 
428 (1971). 

[12] We conjecture that the ferromagnetic fluctuation is the 
most advantageous of spin ordered OFs. It gives a 
higher DOS than the paramagnetic OF, see I. Vekhter, 
A. V. Shytov, I. A. Gruzberg, and A. V. Balatsky, 
cond-mat/0210421 to appear in Physica B. 



